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Abstract. 

We study effects of electron-electron interactions and confinement potential 
on the magneto-optical absorption spectrum in the far-infrared range of lateral 
quantum dot molecules. We calculate far-infrared (FIR) spectra for three 
different quantum dot molecule confinement potentials. We use accurate exact 
diagonalization technique for two interacting electrons and calculate dipole- 
transitions between two-body levels with perturbation theory. We conclude that 
the two-electron FIR spectra directly reflect the symmetry of the confinement 
potential and interactions cause only small shifts in the spectra. These predictions 
could be tested in experiments with nonparabolic quantum dots by changing 
the number of confined electrons. We also calculate FIR spectra for up to six 
noninteracting electrons and observe some additional features in the spectrum. 
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1. Introduction 

Far-infrared (FIR) magneto-optical absorption spectroscopy is one experimental 
technique to study electrons confined in semiconductor quantum dots (QDs) [3,0. 
It was, however, realized at the early stage of QD research that FIR spectroscopy is 
unable to reveal interesting many-electron effects in parabolic-confinement QDs. This 
is because the electromagnetic waves couple only to the center-of-mass variables of 
electrons. The resulting FIR spectrum is rather simple, showing two branches, o;±, as 
a function of magnetic field [slti • These branches, one with positive energy dispersion 
and one with negative energy dispersion as a function of magnetic field, are 
called the Kohn modes. The spectrum does not depend on the number of electrons 
nor on the interactions between them. This condition in parabolic QDs is called the 
generalized Kohn theorem 5]. The condition can be lifted if the confinement is not 
parabolic, and many experiments show more complex FIR spectra [l|, |y, L1l|. Also 
calculations [HI El [13, El E3 E3 have shown non-trivial FIR spectra of non- 
parabolic QDs. It is also possible that spin-orbit interaction El an d impurities near 
quantum dots 18] can have an effect on the FIR spectrum. In a nonparabolic QD, the 
relative internal motion of electrons could be accessible with the FIR spectroscopy, but 
recent studies suggest that additional features in the FIR spectra are still of collective 
nature [2j- The interpretation of the observed FIR spectra of a nonparabolic QD is 
usually far from trivial. It is clear that deviations arise from nonparabolic confinement, 
but the detailed cause of the deviations, thus the interpretation of spectrum, is 
not always straightforward. It is especially interesting to see how many-electron 
interactions appear in the FIR spectra of QDs. 

In the calculation of FIR spectra of two-electron lateral double QD [3 and lateral 
four-minima quantum dot molecule (QDM) yJJ it was shown that the deviations from 
the Kohn modes arise mainly from nonparabolic confinement and interactions have 
only minor effects on the spectra. Effects of relative motion of electrons on FIR spectra 
were studied by turning electron-electron interactions on and off. The results indicate 
that the FIR spectra of QDMs reveal mainly the center-of-mass collective excitations 
of electrons. Actually more pronounced deviations from Kohn modes are observed 
when interactions between the electrons are turned off. This is an interesting result, 
suggesting that less features are observed with more electrons in a non-parabolic QD. 
These studies, however, include only two electrons and the FIR spectra may depend 
on particle number. Also as more single-particle levels are occupied the center-of-mass 
excitations with slightly different energies may be observed where interactions do not, 
necessarily, have any significant role. Calculations where electron interactions can be 
turned on and off could show how the electron-electron interactions show up in FIR 
spectra for greater electron numbers (N > 2). 

In this article we study three different lateral quantum dot molecule (QDM) 
confinements. We study two- minima QDM (double dot), and square-symmetric and 
rectangular-symmetric four-minima QDMs. We calculate the FIR spectra for all 
the QDM confinements and compare them to the parabolic-confinement QD FIR 
spectrum. We analyze in detail the effect of electron-electron interactions on the FIR 
spectra. We use very accurate exact diagonalization technique for interacting electrons 
but limit our studies to two interacting electrons. We also calculate the FIR spectra 
of non-interacting electrons up to six electrons by occupying single-particle levels and 
calculating dipole- transit ion elements between the single-particle levels. These provide 
some insight how FIR spectra are modified with greater electron numbers, even if we 
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do not include electron-electron interactions. Especially if excitations are collective 
also with the greater electron numbers, the FIR spectra are basically given by the 
transitions between single-particle levels of the quantum dot confinement in question. 

2. Model and method 

We model the interacting two-electron quantum-dot molecule with a two-dimensional 
Hamiltonian 

i=l x ' 

where V c is the external confinement potential of the QDM taken to be 

V c (r) = l -m*iolmm j^(r - L,)'} . (2) 

The coordinates are in two dimensions r = (x,y), the L/s (Lj = (L x ,L y )) give 
the positions of the minima of QDM potential, and M is the number of minima. 
When Li = (0,0) (and M = 1) we have a single parabolic QD. With M = 2 and 
Li 5 2 = (±1^,0) we get a double-dot potential. We also study four-minima QDM 
(M = 4) with minima at four possibilities of (±L x ,±L y ) (see Fig. [I}. We study 
both square- symmetric (L x = L y ) and rectangular- symmetric (L x ^ L y ) four- minima 
QDMs. The confinement potential can also be written using the absolute values of x 
and y coordinates as 

V c (x,y) = ^m*^ x 

[r 2 -2L x \x\-2L y \y\+L 2 x + L 2 y ] . (3) 

For non-zero L x and L y , the perturbation to the parabolic potential comes from the 
linear terms of L x or L y containing also the absolute value of the associated coordinate. 

We use the GaAs material parameters m*/m e = 0.067 and e = 12.4, and 
confinement strength Tiljq = 3.0 meV. This confinement corresponds to harmonic 
oscillator strength of Iq = ^/h/uom* ~ 20 nm. We concentrate on closely coupled 
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Figure 2. Relative transition probabilities between Fock-Darwin energy levels 
P=(cf> n >, l±1 \e±^r\cf> n j). 



QDMs where L XiV < lo. The magnetic field (in z direction) is included in the 
symmetric gauge by vector potential A. The Hamiltonian of Eq. is spin-free, 
and the Zeeman energy can included in the total energy afterwards (Ez = g* HbBSz 
with g* = —0.44 for GaAs). We disregard the threefold splitting of each triplet state 
(Sz = 0, ±1) and consider only the lowest energy one (Sz = 1)- 

We exclude the explicit spin-part of the wave function and expand the many- 
body wave function in symmetric functions for the spin-singlet state (S = 0) and 
anti-symmetric functions for the spin-triplet state (S = 1). 

^Mri,r 2 ) = ^OLi,j{<j>i{v{)<j)j{T2) 

+ (-l) s 0,(r 2 )0 J (r 1 )}, (4) 

where o^j's are complex coefficients. The one-body basis functions <fii(r) are 2D 
Gaussians. 



(r) = x nx y ny e 



-r 2 /2 



(5) 



where n x and n y are positive integers. The complex coefficient vector a^j of ^\ 
and the corresponding energy E\ are found from the generalized eigenvalue problem 
where the overlap (S) and Hamiltonian (H) matrix elements are calculated analytically 
(Hai = EiSai). The matrix is diagonalized numerically. 

The FIR spectra are calculated as transition probabilities from ground state 
(£"o) to excited states (E{) using the Fermi golden rule within the electric-dipole 
approximation: 

\ 2 



Ai± oc 



9i 



5{Ei -E - M- 



(6) 



We assume circular polarization of the electromagnetic field: e ±l ^^ i Y i = =^ 
iyi) = z±, where plus indicates right-handed polarization and minus left-handed 
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Figure 3. Transitions for spin singlet (a) and spin triplet (b). (c) Two-particle 
energy levels of a parabolic quantum dot. Red lines show singlet energy levels 
and blue lines triplet energy levels. Gray background shows singlet ground state 
regions and white triplet ground state regions. Dipole transitions from one spin 
type to another are forbidden. 



polarization. The results are presented for non-polarized light as an average of the 
two circular polarizations. 

3. Far-infrared spectra of two-electron quantum-dot molecules 

The dipole-allowed magneto-optical excitation spectrum of an isolated harmonic- 
confined QD consist of two branches cj + and c^_, whose energy dispersion is well 
understood and does not depend on the number of electrons in the quantum dot: 

AE± = huo± = fi^Juol + (uo c /2) 2 ± huo c /2. (7) 

Mo describes the external confinement, uo c = eB/m* is the cyclotron frequency, and m* 
is the effective mass of electron. Fock-Darwin energy levels as a function of magnetic 
field are given by E n \ — (2n+ |Z| + l)Tiuo — \lTiuo c , where n = 0, 1, 2, ... and I = 0, ±1, ... 
are principal and azimuthal quantum numbers, respectively, and uo = y / 'uj^ + (co c /2) 2 . 
In the dipole-allowed transitions between Fock-Darwin states angular momentum must 
change by unity, Al = ±1. Fig. shows relative transition probabilities between Fock- 
Darwin energy levels P = (0 n ',z±i |e ± ^r|^ n ,/). 

For a two-electron QD the ground state can be either spin singlet (S — 0) or spin 
triplet (S = 1) depending on the magnetic field strength. Even if many-body energy 
levels have more complicated magnetic field dispersion than the single-particle levels, 
the dipole-transitions in the many-body case always equal to transitions between the 
single-particle levels in a parabolic confinement. A schematic picture of transitions at 
zero magnetic field is shown in Fig. 0(a) and (b) for spin-singlet and spin-triplet states, 
respectively. Fig. (c) shows a few lowest two-body energy levels as a function of 
magnetic field for the singlet and triplet states of a parabolic QD. The gray background 
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marks the magnetic field region of the spin-singlet ground state and white the spin- 
triplet ground state. Dipole transitions from one spin type to another are forbidden. 

3.1. Interacting two- electron spectra 

Fig. 21 shows the calculated far- infrared (FIR) absorption spectra for two interacting 
electrons. Fig. 0] ( a ) shows the FIR spectrum for parabolic QD, (b) for two- 
minima QDM (double dot), (c) for square symmetric four- minima QDM and (d) 
for rectangular- symmetric four-minima QDM. The energy of absorbed light is given 
in meV as a function of magnetic field. The width of the line is proportional to the 
transition probability, which is also plotted below each spectrum (in arbitrary units) 
for each of the branches in the spectrum. There are, of course, many more possible 
transitions with nearly zero or small probability. Only the lines with transitions 
probability exceeding 1% of the maximum transition probability are included in Fig.QJ 
Large open circles in the spectra represent the two Kohn modes, uo + and uo- , plotted 
with IT spacings. In QDMs there are ground state transitions from spin-singlet 
(5 = 0) to spin-triplet (S = 1) ground states, and the other way round, as a function 
of magnetic field [l9|,|20(. Vertical lines indicate the singlet-triplet (or triplet-singlet) 
transition points and the singlet ground state regions of spectrum are marked with 
red color and the blue lines denote triplet regions. Dipole transitions from one spin 
type to another are forbidden. 

Fig. IH (a) shows the calculated FIR spectra for a single parabolic QD. The dipole 
transition probabilities are calculated from the two-body ground state level to higher 
two-body energy levels as shown in Fig. 01(c). The spectrum of Fig. El (a) shows two 
Kohn modes o;+ and uj-. These coincide perfectly with the open circles presenting 
Kohn modes, as they should, since in a parabolic QD the FIR spectra does not depend 
on the number of electrons in the QD nor on the interactions between them p. The 
magnitude of the absorbed light can only change when the number of electrons in the 
QD is changed |l|, |4( . However, in the two-electron QD of Fig.0](a) the width of the 
upper branch is halved after the first singlet triplet transition. The transition line is 
split into two degenerate transitions with almost equal transition probabilities, which, 
if summed up, would equal to the transition probability of uo-. In the calculations we 
can identify two degenerate levels to which electrons can be excited (see Fig. 01(b)), 
but in experiments the observed quantity is the overall absorption of the energy in 
question. Therefore one can observe Kohn modes with constant transition probabilities 
as a function of magnetic field. 

Fig. EI (b), (c), and (d) show FIR spectra for quantum-dot molecules. Now as 
the symmetry of the confinement is lower, the center of mass and relative motion do 
not decouple. There are clear deviations from the Kohn modes. In all of the QDM 
potentials we can see some general features. The main branches of the spectra always 
lie below the Kohn modes (below open circles). There are clear anticrossings and the 
upper branch is split to two in some parts of the triplet spectra. 

In the double dot spectra of Fig. 01(b), the excitations at zero magnetic field are 
identified to a lower energy one along the long axis (the line connecting two dots, x 
axis) and a higher energy one along the short axis (y axis) . This can be also verified by 
calculating the absorption of linearly polarized light 9] . At B = the excitation of u+ 
is unaffected by the interactions between electrons and coincides with the confinement 
{Tiuoq = 3 meV), because in the double dot the potential is still parabolic along the 
y axis. However, with non-zero B this is no longer true since a magnetic field mixes 
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Figure 4. Calculated far-infrared spectra of single parabolic quantum dot in (a), 
double dot with L x = 15 nm in (b), square-symmetric four-minima quantum- 
dot molecule with L x = L y = 10 nm in (c), and rectangular-symmetric four- 
minima quantum-dot molecule with L x = 5,L y = 10 nm in (d). The energy 
of absorbed light is given in meV as a function of magnetic field. The width of 
the line is proportional to the transition probability, which is also plotted below 
each spectrum (in arbitrary units) for each of the branches in the spectrum. 
Transitions with small probability are excluded from the spectra. Only the lines 
with transitions probability exceeding one percent of the maximum transition 
probability are included in the spectra. Large open circles in the spectra represent 
two Kohn modes plotted with one tesla spacings. Vertical lines indicate the 
singlet-triplet (or triplet-singlet) transition points. 
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the two linear polarizations. At high magnetic field the spectrum approaches Kohn 
modes and no anticrossings are observed. In the high field region the electrons are 
effectively more localized into individual dots [2(j. Then the non-parabolic nature of 
the confinement potential has a less important role, since the electron density is more 
localized close to the parabolic minima. 

After the singlet-triplet transition there appears an additional level, cj + 2, above 
U-\- with a lower transition probability. A simple explanation for two positive dispersion 
levels would be in the picture of Fig. 01(b) where the two transitions of ou+ mode could 
have slightly different energy difference in a non-parabolic confinement. However, the 
interpretation is not so straightforward, since the dipole transitions are calculated 
between two-body levels and not occupied single-particle levels, as we will see when 
we analyze interactions in the next subsection. Another deviation from the Kohn 
modes of a parabolic QD are seen as small discontinuities in the energy of the spectral 
lines when the singlet changes to triplet. However, this discontinuity is very small and 
may not be resolved in experiments. See Ref. 9] for more details on double dot FIR 
spectra. 

In the square-symmetric four-minima QDM (Fig.0](c)) the potential is identical 
in x and y directions resulting in only one excitation at B = as the long and short 
axes are equal. This excitation is clearly lower than the Tujq = 3 meV confinement 
energy of a single minimum. The upper branch is split to two levels in the first triplet 
region at magnetic field values between 1 and 3 T, but in the second triplet region 
B > 4 T there is only one upper branch. In square-symmetric four-minima QDM 
the two levels of the split-off upper branch have almost equal transition probabilities. 
Discontinuities in the spectra appear in the singlet-triplet (or triplet-singlet) transition 
points. A brief representation of the FIR spectra of square-symmetric four-minima 
QDM can also be found in Ref. 

The spectrum of rectangular-symmetric four- minima QDM, in Fig. 0] (d), has 
similar features as the double dot spectra. There is a gap between the two branches at 
B = 0. The lower branch at B = in Fig. 01(d) corresponds to excitation along y axis 
and the upper branch along x axis. The upper mode does not start at 3 meV, as in the 
double dot, since the potential is not parabolic in x direction. After the singlet-triplet 
transition an additional mode, cj + 2, appears above lu+. The transition probability of 
cj + 2 is higher compared to the cj + 2 of double dot, but the transition probability is not 
as high as in the square-symmetric four-minima QDM. The modes of Fig. 0] (d) at 
high magnetic field do not become so close to Kohn modes as in the double dot. The 
reason is that electrons localize into two decoupled double dots rather than into single 
QDs as in the double dot spectrum of Fig.0](b). 

All QDM confinements show anticrossings in uj+. The clearest anticrossings are 
seen at low magnetic field strengths. At higher B there are still anticrossings but 
the energy gap is so small that they are not visible in the spectra. The transition 
probabilities in the lower panels, however, reveal anticrossings as the transition 
probability of one mode decreases and the other one increases as a function of 
magnetic field. The biggest anticrossing gap in energy is seen with square-symmetric 
four-minima QDM. This is interesting, since in our previous studies the square- 
symmetric four-minima QDM resembled the most a parabolic QD whereas double 
dot and rectangular four-minima QDM showed greater deviations in the ground state 
properties and in the low- lying eigenstates |2£|. Obviously this is not true with the 
excitation spectrum and with the higher eigenstates. 
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Figure 5. Far-infrared spectra of two interacting electrons in the singlet 
symmetry in (a)-(c) for double dot, square-symmetric four- minima quantum- 
dot molecule, and rectangular-symmetric four-minima quantum-dot molecule, 
respectively, (d)-(f) show corresponding non-interacting far-infrared spectra. 



3.2. Analysis of interactions in two-electron FIR spectra 

We conclude that different types of deviations appear in the spectrum when the 
confinement potential is not perfectly parabolic. In order to analyze these deviations 
in more detail, we plot the full spectra up to 6 T magnetic field for both singlet and 
triplet states with two interacting electrons and also for two noninteracting electrons 
for the comparison. We are interested to see which features in the spectra are due 
to electron-electron interactions and which result from the lower symmetry of the 
confinement potential. It is also interesting to see if some features in the spectra can 
be used to identify the confinement potential. This would be useful for experiments 
where the profile of the underlying confinement potential is not clear. 

We will first analyze the two-electron singlet spectra shown in Fig. [51 (a)-(c) 
for interacting electrons and in (d)-(f) for noninteracting electrons. The two-electron 
triplet spectra are shown in Fig.EJ(a)-(c) and the corresponding noninteracting spectra 
in (d)-(f). In a real two-electron QDM system, singlet-triplet transitions must be 
taken into account and therefore experimentally observable spectra for two interacting 
electrons are shown in Fig.0] As the non-interacting singlet spectra, in Fig. 03 (d)-(f), 
are the same as the single-particle spectra, these would correspond to experimental 
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FIR spectra with one electron in the QDM. The non-interacting triplet, on the other 
hand, corresponds the FIR spectra of two occupied single-particle levels. 

The double dot singlet spectrum in Fig. [51 (a) looks qualitatively very similar to 
the noninteracting spectrum in Fig. 0(d). The first difference is that the Kohn modes 
lie at the lower energy in the non-interacting spectra compared to the interacting 
spectra. The only exception is in the upper mode at B = where both have excitation 
energy of Tiuoq = 3 meV. In zero magnetic field this corresponds to linearly polarized 
excitation along y axis where the potential is parabolic and therefore excitations 
are not affected by the electron-electron interactions 9]. Another difference is the 
zero-field gap between Kohn modes which is clearly greater in the non-interacting 
case. Higher energy excitations in the interacting case can be explained with the 
Coulomb repulsion between the electrons. Coulomb repulsion effectively steepens the 
confinement resulting in higher excitation energies. The anticrossing points are also 
affected by interactions whereas the gap remains almost the same. The anticrossings 
can be seen at lower magnetic field values in the interacting spectrum. Qualitatively 
the two spectra look very similar and we conclude that the deviations observed in the 
spectrum result from the low symmetry confinement and interactions only shift the 
excitation energies and change the anticrossing points. 

Generally the same conclusions drawn for the double dot hold for the square- 
symmetric four-minima QDM singlet spectra in Fig. (b) and (e). Now in the 
square-symmetric four-minima QDM the zero field excitation is degenerate as the two 
perpendicular directions (x and y) have identical confinement profiles. Interactions 
shift the zero field excitation to higher energy due to Coulomb repulsion. The first 
anticrossing is again seen at lower B in the interacting spectrum. However, this time 
the anticrossing gap is somewhat greater in the interacting case. As a new feature, 
compared to the double dot, there are now discontinuities in the interacting spectrum. 
The discontinuities in (b) are at the transition points when two singlet states with a 
different symmetry cross. At B « 1.8 T the lowest-energy singlet state changes from 
no- vortex state to two- vortex state and the discontinuity at B « 5.2 T is at the 
crossing point of two- vortex and four- vortex states (see Ref. 20] for more details on 
vortex states). However, the discontinuities are seen at the magnetic field values where 
the true ground state is triplet and therefore the discontinuities in Fig. (b) would 
not be observed in FIR experiments. 

In the rectangular-symmetric four- minima QDM (Fig. ( c ) an d (f)) the 
confinement profile is non-identical and non-parabolic in both perpendicular 
directions. Therefore both zero-field excitations are lower than 3 meV and there 
is also a gap between the two modes. Again excitation energies as a function of 
magnetic field are higher in the interacting spectrum and anticrossings are seen at 
the lower magnetic field values compared to the noninteracting spectrum. Interesting 
point, when comparing interacting and noninteracting spectra, is that the zero-field 
excitation energy of the lower mode shifts to 0.30 meV higher and the higher mode 
shifts only 0.09 meV higher when the interactions are turned on. It means that the 
Coulomb repulsion is over three times more significant in the excitation along the long 
axis (y axis) than it is along the short axis (x axis). Electrons are localized into two 
distant double dots. Coulomb repulsion is more important between the two double 
dots as it is inside one double dot, i.e., along long and short axes. 

Next, we do similar comparisons for the triplet spectra of QDMs. The interacting 
spectra are shown in Fig. (a)-(c) and non-interacting in Fig. (d)-(f). We will 
first compare double dot interacting triplet spectra of Fig. (a) to the noninteracting 
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Figure 6. Far-infrared spectra of two interacting electrons in the triplet 
symmetry in (a)-(c) for double dot, square-symmetric four- minima quantum- 
dot molecule, and rectangular-symmetric four-minima quantum-dot molecule, 
respectively, (d)-(f) show the corresponding non-interacting far-infrared spectra. 



spectra of Fig. ED(cl). In general, the same analysis applies to triplet spectra as for 
the singlet spectra of double dot. However, in the triplet state there is an additional 
mode, cj + 2, above the main branch, in Fig.^Ja). The additional mode has clearly 
weaker transition probability than the cu+ mode. In the noninteracting case the upper 
mode is split to two modes which both have equal transition probabilities. In the 
interacting spectrum the upper mode has lower transition probability and soon after 
4 T it completely vanishes from the spectrum. 

In the square-symmetric four-minima QDM also the interacting triplet spectrum, 
as well as the non-interacting spectrum, has two upper modes with almost equal 
transition probabilities (Fig. 0(b)). In four- minima QDM the split-off branch suddenly 
changes to one mode at B « 3.5 T. There the triplet ground state changes from the 
one- vortex to three- vortex solution (see Ref. |20| for more details on vortex states). 
Now the two triplet states actually cross and there occurs a discontinuity in the FIR 
spectrum. Yet as the singlet is the true ground state of the square-symmetric four- 
minima QDM between B w 3 and 4 T, the discontinuity of Fig. (b) cannot be 
seen in experiments. In triplet spectra there are much smaller differences between the 
excitation energies of interacting and noninteracting electrons compared to the singlet 
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spectra. In the singlet state electrons of opposite spins can occupy the same single- 
particle levels in the many-body configurations and therefore Coulomb repulsion is 
more significant for the singlet state. In other words, the Pauli exclusion principle 
keeps the same-spin electrons further apart and therefore Coulomb repulsion has a 
lesser effect on the triplet state compared to the singlet state. 

In the rectangular-symmetric four-minima QDM, the zero-field gap is really small 
and not visible in the interacting spectra of Fig. El (c) and neither in the non-interacting 
triplet spectra of Fig. 0(f). The transition probabilities of the two main branches 
and 00— ) , in the interacting spectrum, are almost equal at zero field, but soon after 
the magnetic field strength increases the transition probability of ou+ decreases at the 
same time as the transition probability of uo+2 increases. After B = 4 T the starts 
to weaken again as the symmetry of the triplet ground state is changing. The ou+2 
dies out continuously. The two upper modes of the non-interacting spectrum have 
almost equal transition probabilities. The interacting modes are higher in energy, also 
at B = 0, and anticrossings shift to lower B in the interacting spectrum. Anticrossing 
gaps are not noticeably affected by the interactions. 

4. FIR spectra of noninteracting electrons 

A general conclusion drawn above is that the noninteracting spectra show more 
deviations from the Kohn modes than the interacting spectra. All the same features 
are present in noninteracting and interacting spectra. In addition, the non-interacting 
spectra show in all cases similar structure as the interacting spectra. Interactions only 
shift some of the features seen in the spectra. Therefore we calculate FIR spectra of 
non-interacting electrons up to N = 6 electrons in the QDMs to see if more electrons, 
even if non-interacting, produce changes in the FIR spectra. 

Before presenting the FIR spectra we plot single-particle energy levels of all three 
QDM confinements in Fig. [3 The solid lines show QDM energy levels and the dotted 
lines Fock-Darwin energy levels of a parabolic QD. In QDMs the energy levels shift 
to lower energies compared to parabolic QD. Also many anticrossings and zero-field 
splittings of energy levels are visible in QDMs. 

We obtain noninteracting FIR spectra of QDMs by occupying lowest single- 
particle energy levels with N non-interacting electrons and calculating dipole 
transitions to higher unoccupied levels. The noninteracting FIR spectra are shown 
in Figs. [SI OH and El for a double dot, square-symmetric four- minima QDM and 
rectangular-symmetric four-minima QDM, respectively, (a) to (d) correspond to FIR 
spectra of N = 1 to N = 6 electrons confined in a QDM. N = 1 and N = 2 are 
also shown in Figs. [51 and [SJ The transition probabilities are multiplied when more 
electrons are added in the QDM. We divide each transition probability with N in 
order to ease comparisons. 

We will first analyze the single-particle spectra of a double dot. Fig. [HI (c) for 
N = 3 shows an additional mode below uo- in the low-field region with the excitation 
energy below E — 2 meV. Also another mode is visible with a small transition 
probability which has zero excitation energy at B « 0.9 T. If we study the single- 
particle energy levels of Fig. Q (a), one can see that with the three lowest levels 
occupied, the uppermost levels cross at B « 0.9 T. When these levels cross, the weak 
mode has excitation energy of AE = and the transition probability with the E < 2 
meV mode vanishes. Another feature that is observed with N = 3, but absent in 
two-particle spectra, is an anticrossing in cj_. This anticrossing is opposite to the 
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Figure 7. Single-particle energy levels of double dot (L x = 15, L y = nm), 
four-minima square-symmetric quantum-dot molecule (L x = L y = 10 nm) and 
four-minima rectangular symmetric quantum-dot molecule (L x = 5, L y = 10 
nm). Dotted lines show single-particle (Fock-Darwin) energy levels of parabolic 
(hujQ = 3 meV) quantum dot. 



anticrossings in uo+ where the low-field mode curves upwards and the high-field mode 
curves from down to up while increasing its strength. cj + with three electrons looks 
quite similar as 00+ in N = 2 spectrum. 

In N — 4 spectrum of Fig. |HI (d) there is one level below uo- at low B but the 
clear anticrossing, present in N = 3, is missing from the spectrum. After B w 1.7 
T there is just one lower mode. The N = 5 and 6 spectra are rather featureless as 
clear anticrossings and additional modes are missing. There are many different levels 
but as they are lying close in energy, the overall spectrum resembles just two Kohn 
modes without clear additional features. In general, anticrossings are smaller and fewer 
clearly discernible additional modes are visible with N = 5 and 6. Actually, there is 
some structure of uo- of N = 5 and N = 6 but as energies lie so close, the structure 
would be very difficult to observe experimentally. However, if the perturbation from 
parabolic potential is stronger these levels may have larger energy separation as in the 
experiments of Hochgrafe et al p| . Another noticeable feature in the noninteracting 
FIR spectra of double dot is the zero field gap between the two main modes. For 
N = 1 it clearly has the greatest value whereas for other electron numbers it does not 
vary much. 

Fig.Elshows non-interacting FIR spectra of square-symmetric four-minima QDM. 
In N = 3 spectrum of Fig. Efc) there is one additional mode below uo- and one 
anticrossing in uo- . The additional mode below uo- is much closer to uo- in energy as 
it is in the double dot spectrum. Also u- has lower transition probability than the 
additional mode in four-minima QDM. N = 4, 5 and N = 6 spectra of four-minima 
QDM show multiple modes of cj + and some structure in cj_ , but these are very close 
in energy and may not be resolved in experiments. However, there is a discernible 
split-off mode below 00+ mode with a lower transition probability. It is difficult to say 
how anticrossing gaps of cj + at B w 1.7 T are changing with increasing N because of 
multiple modes. One interesting thing is that zero-field excitation of the degenerate 
main branches is higher, E w 2.6 meV, with N = 3 — 6 than it is with N = 1 and 2, 
where E ^ 2.1 meV. 




Noninteracting FIR spectra of rectangular-symmetric four-minima QDM in 
Fig. El does not deviate significantly from the FIR spectra of two other QDM 
confinements. Additional mode below uo- is very weak but visible and exist only 
at very low magnetic field strengths in rectangular-symmetric four-minima QDM for 
TV = 3, 4 and N = 6 in Fig. EH Clear anticrossing in cj_ occurs only with N = 3. 
The zero-field gap between cu+ and u- becomes very small with with more than one 
electron in the QDM. After the first anticrossing, at B « 2 T, the uj-\. is split to two 
modes with N > 1, where the lower mode has lower transition probability. 

5. Summary and conclusions 

In the calculated far-infrared (FIR) spectra we are able to compare three different 
quantum-dot molecule (QDM) confinement potentials and study the effects of 
interactions on the two-electron FIR spectra. We have also calculated FIR spectra 
for a few noninteracting electrons. The three studied lateral QDM confinements are 
a two-minima QDM (double dot), and square-symmetric and rectangular-symmetric 
four-minima QDMs. 

We observe anticrossings in the upper cj + mode for both interacting and 
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Figure 9. Far-infrared spectra for N = 1 — 6 non-interacting electrons in square- 
symmetric four-minima quantum-dot molecule (L x = L y = 10 nm). 



noninteracting spectra. The anticrossing gaps are not significantly altered by the 
interactions, but the anticrossing positions are shifted to lower magnetic field values 
with the interacting electrons. When we increase the particle number in the 
noninteracting FIR spectra, the anticrossings in modes become less clear as 
the main mode is split to multiple modes at the anticrossing point. With three 
noninteracting electrons (N — 3) we see a clear anticrossing also in the lower uj- 
mode in all three QDM confinement potentials. 

The uj-\- mode is split to two for N > 2 noninteracting electrons and also in the 
interacting spectra when the spin-triplet is the ground state. With two interacting 
electrons this additional mode vanishes from the spectrum at greater magnetic field 
values, approximately at B > 4 for all QDM confinements. The upper split-off mode 
has a lower transition probability in the interacting case in double dot and rectangular 
four- minima QDM confinements. In the noninteracting case, even if we add electrons 
up to six electrons in a QDM, the o; + consists of only two energy-split excitations, 
where one of the modes, typically the lower one, can have smaller transition probability 
in all three QDM confinements. 

The third point that we have analyzed in the FIR spectra is the splitting between 
the main modes, oj±, at B = 0. The splitting is seen for two rectangular confinements, 




double dot and rectangular-symmetric four-minima QDM, where the excitations along 
short and long axes are not equal. In the two rectangular confinements, the greatest 
gap is seen for single electron. For all N > 2 the gap is clearly smaller compared 
to N = 1 but does not change significantly between N = 2 — 6 noninteracting 
electrons in the QDM. The zero-field splitting usually diminishes as electron-electron 
interactions effectively steepen the potential which leads to increase in excitation 
energies. Interactions have stronger effect on the lower mode, which corresponds 
to excitation along the long axis of the confinement, therefore decreasing the gap. 

Generally, QDM confinements have lower excitation energies compared to a 
parabolic quantum dot with the same confinement energy Tiujq. When interactions 
are turned on, the Coulomb repulsion effectively steepen the potential leading to a 
small increase in the excitation energies. 

The same features are in general present in interacting and noninteracting spectra, 
but small shifts may occur with interacting electrons compared to the noninteracting 
FIR spectra. Therefore, based on calculations with two interacting electrons, we 
conclude that deviations from the Kohn modes are due to nonparabolic confinement 
potential. Interactions only shift some features in the observed spectra but otherwise 
excitations are of collective nature. However, this study is only for two interacting 
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electrons and maybe more features can be seen with more interacting electrons in 
a nonpar abolic quantum dot confinement. Also it is possible, but unlikely, that 
different nonparabolic confinement potentials could show more interaction effects in 
FIR spectra. 

We have studied how the FIR spectra change when interactions are turned on 
and off. These predictions could be checked in experiments when the number of 
electrons is increased in the dots, even if the limit of single or two electrons could not 
be achieved. These predictions are that i) anticrossings are shifted to lower magnetic 
field values, ii) excitation energies increase and Hi) zero field gaps decrease as the 
Coulomb interaction increases in the quantum dot (more electrons inside dots). It is 
also possible that the confinement potential may change as the number of electrons 
changes in the dot [2l[. In a rectangular symmetry also polarized light can be used 
to analyze excitations. These features are seen, in general terms, in elliptic quantum 
dots studied by Hochgrafe et al However, they observe a bit more complicated 
FIR spectra with multiple negative and positive dispersion Kohn modes. 

This work has been supported by the Academy of Finland through its Centers of 
Excellence Program (2000-2005). 
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